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The image information acquisition ability of a conventional camera is usually much 
lower than the Shannon Limit since it does not make use of the correlation between 
pixels of image data. Applying a random phase modulator to code the spectral images 
and combining with compressive sensing (CS) theory, a spectral camera based on 
true thermal light ghost imaging via sparsity constraints (CISC spectral camera) is 
proposed and demonstrated experimentally. CISC spectral camera can acquire the 
information at a rate signihcantly below the Nyquist rate, and the resolution of the 
cells in the three-dimensional (3D) spectral images data-cube can be achieved with a 
two-dimensional (2D) detector in a single exposure. For the hrst time, CISC spectral 
camera opens the way of approaching the Shannon Limit determined by Information 
Theory in optical imaging instruments. 


sshan@mail. shcnc .ac.cn 


1 





I. INTRODUCTION 


Conventional Camera, as one of the most important appliances to get image informa¬ 
tion, records the image of an object based on the point-to-point correspondence between 
the object-space and the image-space. Becanse the correlation between pixels of image^cant 
be applied, the image information acqnisition efficiency of snch conventional point-to-point 
imaging mode is much lower than the Shannon Limit^>^ determined by Information Theory 
in optical imaging instruments^^—. Unlike the conventional direct point-to-point imaging 
mode, the resolution of the pixels of ghost imaging is determined by the correlation of light 
held huctuations corresponding to the two pixels respectively, which can be measured on-line 
or pre-determined.— Combining with compressive sensing (CS) theoryi^i^'— , ghost imaging 
via sparsity constraints (CISC) has many potential applications including super-resolution 
imagingl^^— , three-dimensional (3D) computational imaging with single-pixel detectors^^, 
3D remote sensing^i^i^ imaging through scattering mediai^^i^^, object tracking^i, object 
authenticatiort^^i^ and X-ray Fourier transform diffraction imaging^Si^— . For thermal light 
ghost imaging, according to the illumination source, it can be classihed to two categories: 
ghost imaging with pseudo-thermal light and true thermal light. Ghost imaging with true 
thermal light and sunlight have been respectively demonstrated by detecting the temporal 
fluctuation of thermal light and applying the intensity correlation between the intensity dis¬ 
tributions at the reference arm and the test arm.— Comparing with ghost imaging with 
pseudo-thermal light, this scheme of ghost imaging with true thermal light has to face the 
difficulty of detecting the temporal fluctuation of true thermal light which requires the re- 
spouse time of detector less than the coherence time of true thermal light r = —— oc ^ 
( A is the wavelength, AA is the linewidth, c is the speed of light) which can be as short as 
femtosecond. In order to increase the coherence time of the illumining true thermal light, 
monochrome imaging is required which results in the vast majority of radiation energy from 
the target scene being hltered out, making the energy efficiency of ghost imaging applying 
the temporal fluctuation of true thermal light very low. Moreover, the fluctuating true ther¬ 
mal held needs to be split before the light held illuminating the object in the system and 
recorded in the reference pathwhich makes the scheme even more difficult to be applied in 
remote sensing. In this paper, for the hrst time, we propose a spectral camera based on 
true thermal light ghost imaging via sparsity constraints (CISC spectral camera) without a 
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splitter. CISC spectral camera modulates the true thermal light into a spatially fluctuating 
pseudo-thermal light using a spatial random phase modulator— >21 which, at the same time, 
also acts as a random grating generating the uncorrelated speckles for different wavelengths, 
the 3D spectral images data-cube is then modulated into a two-dimensional (2D) data plane 
and CISC spectral camera can achieve the whole wavelength image in a single exposure, 
leading to a more convenient detection process and higher energy efficiency compared to 
ghost imaging applying the temporal fluctuation of true thermal light. Combining with CS, 
CISC spectral camera can acquire the information at a rate significantly below the Nyquist 
rate which opens the way of approaching the Shannon Limit determined by Information 
Theory in optical imaging instruments^'-'^. 


II. SCHEMATIC & RESOLUTION 


The schematic of CISC spectral camera is shown in FigUl The system consists of (1) 
an imaging system, which projects the object image in the object plane ‘a’ onto the first 
image plane ‘6’, (2) a spatial random phase modulator, which disperses the image with 
different wavelengths as a random grating and modulates the image to generate the speckles 
in plane ‘c— (3) a microscope objective, which magnifies the speckles in plane ‘c’, and 
(4) a charge-coupled device (CCD) detector recording the magnified speckles. 
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FIG. 1. Schematic of CISC spectral camera, (a) The object plane; (b) the first image plane; (c) 
the speckles plane; (1) an imaging system; (2) a spatial random phase modulator; (3) a microscope 
objective; (4) CCD detector. 


Denoting the spectral light intensity distribution in the first image plane ‘b’ by A (r*. A;) 
and the intensity distribution in plane ‘c’ by U (rt) respectively, we have^ 


Ur,)=//k{r„X,)h,(rur„X,)dr.dX,. 


( 1 ) 
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where hi{rt]ri, Xi) is the incoherent intensity impulse response function, rt is the coordi¬ 
nate in plane ‘c’, r* and A; are respectively the coordinate and wavelength of the light 
intensity distribution in the first image plane ‘6’. To record the pre-determined refer¬ 
ence spatial intensity fluctuation of the pseudo-thermal light without objects, a coherent 
monochrome point source at pixel r' with wavelength A; in the first image plane ‘6’, de¬ 
noted as hri'f'i, M', X'l) = 5{ri — r'. A; — A;), is used to illuminate the spatial random phase 
modulator, and the recorded light intensity K) f^e plane ‘c’ is given by 

hr{ruri,X'i) = J j Ib^{ri,Xi)hj{rt;ri,Xi)dridXi 

= j J 6{ri-rlXi-X'i)hi{rt;ri,Xi)dridXi ( 2 ) 

= hi{rt-,rl,X'i). 

During the imaging process, the intensity distribution in the hrst image plane ‘6’ ^i) 

is simply the image, denoted as Ti{ri, A;), of the object Tg^rg, Xi) in the object plane ‘a’. 


htir^Xi) = Ti{ri,Xi). 


(3) 


Combining Eqs. ([I])© with ([3]), the intensity distribution in the speckle plane ‘c’ is 

^ctin) = J J Ti{ri,Xi)Ic,{rt;ri,Xi)dridXi. ( 4 ) 

Eq. (jl]) shows that is the Ti{ri, Xi) weighted integration of the pre-determined reference 

spatial intensity fluctuation of pseudo-thermal light Xi). Therefore, each pixel 

of CCD detector is equivalent to a measurement of the bucket detector in the test arm of 
ghost imaging scheme. The second-order correlation function between the spatial intensity 
fluctuation in the pre-determined reference arm and test arm is dehned as 


a;) = {Eliru r', A;)B;(r,)B„(r,)B.„(r,;r', A;)).., 


( 8 ) 


where is the ensemble average about the coordinate of the light intensity distribution 

Tt- Combining Eqs. ©(jl]) with (|5|), the second-order correlation function Aj) is given 

by 


G‘"’(>-'.A;) = I |r,(r.,A,)Gg>(r„A,,r',Ai)dr.A. 


( 6 ) 


where A^, r', X\) = r', X'i)E*^{rt] n, Xi)EcXrt] n, Xi)EcXrt] r', A]))^, is the second- 

order correlation function of the light helds at different pixels and wavelengths in the hrst 
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image plane ‘ 6 ’. In order to calcnlate G® (rj, A/, r', A^, the height antocorrelation fnnction 
of the spatial random phase modulator is assumed as^ 

_ / 2 

Rvi^o,ro) = {v{ro)v{ro)) = ^ = i?^(Aro), A(ro) = ro - r^, (7) 

where //(rg) and ^ 7 (^ 0 ) are respectively the height of the spatial random phase modulator at 
tq and Tq, uj and ( are respectively the height standard deviation and transverse correlation 
length of the spatial random phase modulator. Assuming that the light field fluctuations in 
the speckles plane ‘c’ corresponding to pixel r' in the hrst image plane ‘ 6 ’ with wavelength 
A; obeys the complex circular Gaussian distribution, X'l) can be written as 


where 


Gfjiru A«,r', A^) = X'i)El{rt;ri, Xi)Ec,irt;ri, Xi)EcXrt]r[, X'i))r, 

= {EliruK, X[)E^Xrur[, A;))^,(E*^(rt;r*, Ai)^c,(G;ri, AO)^ 

+ r'i, X'i)El{rt-, n, Az))^,(Ec,(rt; n, Xi)Ec,{rt; r', A^))^, 

“ Xi))rt{GriTti'^ii Xl))rt \ Xl, T^, Xi) \ 

= {Icr{rt-,r-, X[))r,{Icrirt;ri, Xi))r,[l + gfj{{ri, Xi,r[, X[)], 
Jcriri,Xi,r[,X'i) = {El{rt;r[, X'i)El{rt;ri, Xi))r,, 

I JcM,^hr'-,X'i) f 


( 8 ) 


Ai,r', a;) = 


{Icr{rt;r'i,X'i))r,{Icr{rt-,ri,Xi))r,' 


(9) 

( 10 ) 


g^\{ri, Xi,ri, X'l) is defined as the normalized second-order correlation function of the light 
helds at different pixels and wavelengths in the hrst image plane ‘ 6 ’. According to the Fresnel 
diffraction theorem, the light held in the speckles plane ‘c’ propagated from pixel r' in the 
hrst image plane ‘ 6 ’ with wavelength A( is 

exp{j7r[2{zi + Z 2 ) + {n - + z2)]/X'i} 


Ecr{rt]r'i,X'i) = 


jX'iZiZ2 


f WN r.T^{zi + Z2), zin + Z2r'i j 
X / t{rQ,Xi)exp[]—-, -(ro-^-) ]dro. 


( 11 ) 


X'iZiZ2 ' Zi -h Z 2 

where t(ro, A;) = ea;p[j27r(?7, — l)? 7 (ro)/A(] is the transmission function of the spatial random 
phase modulator. (/c,.(rt; r'. A;)) and (/c^(rt; r*, A^)) are respectively given by 
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{IcXrt]r'i,X'i))r, = {ElXrt]r'i,X'i)EcXrt]r'i,X'i))r, = 
(A.(a; A, A;))^, = {E*^{rt]ri,Xi)E^Xn]ri,Xi))r^ = 


A^2A;(^1 + ^ 2 ) 

1 

ZlZ 2 Xi{zi + Z 2 )' 


( 12 ) 

(13) 
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Substituting Eqs. d7|) fl^ (fTTj) into ffTOj) yields 


J2),,r A r' y) = + !)!'(? + ?)) 

9 cr \y ^l) 


( 14 ) 


X I / / exp{[2n{n - ^ - y)}d/idz/p, 

A/A; Z 1 Z 2 A; A; 


where 


/i = ro 


Zir* + 


u = r, 


Zl + X2 

, _ Zin + X2r( 

° Zi + Z2 ' 

AA( = a; - Az, 

u, ,0 ^2(ri-r') 2 t ,,2„_r r/,. ,a ^ 2 Ar '^^^2 


Rriip, v) = (U^expj —[(/i — v) 


Zl + Z2 


-/C] } exp{-[{p-iy) 


Zl + Z2 


Ar' = Vi — r[. 


( 16 ) 

( 16 ) 
(17) 

Kf], (18) 

( 19 ) 


Assuming -C 1, we have 

{zi + Z2yexp{-[27Tu{n - + ]^)} 


ZizlXiX'i 


( 20 ) 


l2 z/) 


X I / / e(rp{[27r(n-l)]2^h-^^}exp{j 


.T^{Zi + Z2), 2 2 , 2 


(/i^ — z/^ H—^z/^)}d/i(iz/p. 
A; 


Assuming S> 1, and the diameter a of the illuminated region in the spatial random 
phase modulator by each cells of 3D data-cube in calibration satishes 7r(j^/A/X2 < 1, 
9 f}{{ri,Xi,r[,X[) is given by 

!lg'((r.,A,,r',A;)«exp{-|27r(*,(n-l)|'^(i-i)nexp{?2^A|;A!{exp{-|^-^5^l2}- 

( 21 ) 

Taking Eqs. fl8|l flT^ flT^ and fl^ into Eq. ([6]), we get the correlation function of intensity 
fluctuations^ 


In 


AG(")(r',A() 

= - {Icr{n]ruXi))rAht{.rt)) 

~ - 727277 ^ 7 ^(a '1 ^i)}< ® {exp{2[2T^uj{n - l)k[f{exp{-[ 
Z{Z^{Zi+Z2Y 

® exp{ — [2Tiu{n — l)]^/i)f}, 


/l 2 r._r r ,2 


{zi + Z2)C 


r}-mr'M 

( 22 ) 


!}}• 
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where fc; = ^, = ^, (g) denotes the operation of convolution. Eq. specihes 

that Ti{rl, X'l) can be separated from the correlation function of intensity fluctuations 
A;), and the resolution is determined by the normalized second-order correla¬ 
tion X'l) at different pixels and wavelengths in the first image plane ‘6’. When 

Vi = r'i, according to Eq. (na, the normalized second-order correlation function of the light 
fields at pixel r- in the hrst image plane ‘6’ with two different wavelengths is given by 


9 fJ{K,Xi,r'i,X[) 

{zi + Z2yexp{-[27TUj{n - 1)]^(4 + ^)} 


.7r{zi + Z 2 )., 
expij --[(ro - 


zlzlXiX'i 

ZiTt + 


X lyyea:p{[27r(n — !)]■ 
ZiTt + Z2r': 


<Rviro,ro) 


Aza; 


} (23) 


2 : 1^2 


--—r/X, - (r' - 'jr /K]}drodr',\\ 
Zi -r Z2 ~r Z2 


Similarly, when Xi = X'l, according to Eq. flT^ . the normalized second-order correlation 
function of the light helds at two different pixels in the hrst image plane ‘6’ with wavelength 
A; is given by 


9fJ{ri,X[,r[,X'i) = exp{-2[27ruj{n - 1)/A;]^{1 - exp[-{ 


Z2^r'- 




(24) 


' {Zi -f Z2)C' 

Fig,2(a) and Fig,2(b) respectively show the comparison of A/, r', A() and g^^{ri^ A(, r', A^ 

between experiment and theory, and the experiment diagram is given in Fig{T] with zi = 
20mm, Z 2 = 0.3mm, u = 2.1pm, ( = 16.75pm, n = 1.516 and the central wavelength 
A; = 600nm. 




(a) (b) 

FIG. 2. (a) The normalized second-order correlation function of the light fields gcj{r[,Xi,r[,X'i) at 
pixel in the first image plane ‘5’ with two different wavelengths; (b) The normalized second-order 
correlation function of the light fields gcj {ri,X'i,r[,X'i) at two different pixels in the first image 
plane T’ with wavelength A^. 
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III. THE MEASUREMENT MATRIX & RECONSTRUCTION 
ALGORITHM. 


There are many methods to improve the imaging quality of ghost imaging.—"— However, 
ghost imaging reconstructions based on the ensemble statistics cannot provide the criterion 
of the necessary number of sampling for a perfect imaging, which makes it impossible to op¬ 
timize the design of ghost imaging system. Combining with CS which provides the recovery 
condition of perfect reconstruction, the quantitative analysis for the necessary measurements 
data can be made. Under the framework of CS theory, the measurement matrix of CISC 
spectral camera is obtained as follows: each of the speckle intensity distributions generated 
by a point light source at pixel r ■ in the A ± AA; spectrum band in the first image plane 
‘6’ is recorded by the randomly selected pixels of CCD detector and reshaped as a 
column vector of length M of the measurement matrix. Repeating the process for all the 
N image pixels in the first image plane ‘6’ and all the L spectral bands, one may have the 
pre-determined random measurement matrix Amxk, where K = L x N. If we denote the 
unknown spectral object image as a iC-dimensional column vector Xkxi, and reshape the 
modulated object intensity distribution recorded by the same M pixels of CCD detector in 
a similar way as a column vector Ymxi, then we may have the discrete from Eq. (HI) . 

Y = AX. (25) 

Spectral object image is usually both spatially and spectrally correlated, which has already 
been utilized in spectral image reconstructions.—"— The reconstruction of the spectral object 
image can generally be regarded as solving a minimization problem which penalizes both 
the Zi norm and the nuclear norm of the data matrix: 

mmniWTpXWi^ + lj.2\\X\\^,s.t.,Y = AX, (26) 

where Xlxn a rnatrix representation of the spectral object image whose columns represent 
different bands of the spectral object image, xjj the sparsifying transform, ni and fi 2 the 
weight coefficients and > 0. In this work, we use a modihed approach based on the 

method described by Eq. 

minpi||?/>X||zj + ij, 2 \\As\\i^, s.t.,Y = AX,X > 0, (27) 
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where As* = Sj+i — si, is the subtraction of the largest singular value si and the other 
Si- The solution of Eq. fl271) tends to have a simultaneous low-effective-rank and sparse 
structure, which much improves the reconstruction quality with low sampling rate. 


IV. EXPERIMENTAL RESULTS.. 

In the experimental setup of CISC spectral camera shown in FigJSl the imaging system 
(Tamron AF70-300mm f/4-5.6) with focal length of / = 180mm projects the object image 
onto the first image plane, a beam splitter (BS) with split ratio 50:50 splits the light field 
into two paths, CCDl detector (AVT Sting F-504C with pixel size of 3.45/rm x 3.45p.m) 
is placed in one of the two paths at the position of the hrst image plane of the system to 
obtain the conventional image of the object for comparison, a spatial random phase mod¬ 
ulator (SIGMA KOKI CO., LTD. DFSQ1-30C02-1000) disperses the images with different 
wavelengths acting as a random grating and modulates the image to generate the speckles, 
a microscope objective with magnihcation = 10 and the numerical aperture N.A. = 0.25 
magnihes the speckles which are then recorded by CCD2 detector (Andor iKon-M) with 
the pixel size 13/rm x 13/rm. The hrst image plane is divided into x Ny = 140 x 140 
pixels with the square of each pixel approximately equal to Ar^ determined by the Eq. fl241) . 
The number of spectrum bands for single exposure is 7, and the images in two wavelength 
ranges of 520 ~ 580nm and 620 ~ 680nm are respectively obtained in two exposures, while 
the theoretical spectral resolution is 20nm in the experimental setup according to Eq. ([23]). 
In order to compare the spectral & spatial resolution of CISC spectral camera with the 
theoretical resolution, as shown in FiglU the spectral object ‘SIOM’ with different parts 
passing through different wavelengths has been selected, and the illuminating source is a 
xenon lamp. The original spectral images of ‘SIOM’ obtained by CCDl detector placed in 
the hrst image plane ‘6’ with corresponding narrowband hlter in front of it are shown in FigH] 
(pixel size is equal to the theoretical resolution of reconstructed images by CISC spectral 
camera for comparing them). The corresponding modulated object intensity distribution Y 
is achieved by CCD2 detector of CISC spectral camera and the reconstructed spectral im¬ 
ages of ‘SIOM’ with 30% sampling rate of 3D date-cube are shown in FigjS] The comparison 
between the original and reconstructed spectral images shows that the resolution of CISC 
spectral camera is in accordance with the theoretical calculation. 
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FIG. 3. Experimental setup of CISC spectral camera. 
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FIG. 4. The original spectral images of ‘SIOM’ obtained by CCDl detector placed on the first 
image plane ‘6’ with corresponding narrowband filter in front of it, showing all the channels from 
620 ~ 680nm. 


The images of the outdoor scene consisting of Mario & Luigi with sunlight illumination are 


shown in Figj6l Fig,6(a) is obtained by a camera, while Fig,6(b) and Fig,6(c) respectively 
show the pictures taken by CCDl detector with narrowband Liters of 550 ± lOnm and 
650 ± lOnm in front of it (pixel size is equal to the theoretical resolution of reconstructed 
images by CISC spectral camera for the sake of comparison). The reconstructed spectral 
images of Mario & Luigi with 30% sampling rate of 3D date-cube are shown in Fig 171 The 
experimental results show that the spectral imaging ability of CISC spectral camera for 
complex scenes is also pretty good. 
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FIG. 5. The reconstructed spectral images of ‘SIOM’ with 30% sampling rate of 3D date-cube, 
showing all the channels from 620 ~ 680nm. 



(a) (b) (c) 


FIG. 6. taken by (a) a camera; (b) GCDI detector passing through narrowband filters of 550 ± 
lOnm; (c) GGDl detector passing through narrowband filters of 650 ± lOnm. 

V. DISCUSSION & CONCLUSION. 

Based on Information Theory, the transmitted information of an imaging system can be 
described by the entropy^*^ 

= - Z!p(u) ^ogp{xi), (28) 

i 

where p{xi) is the probability of Xi occurrence. For the conventional direct point-of-object- 
space to point-of-image-space imaging mode, the conditional entropy H{X\Y) = 0, and thus 
the channel capacity of the conventional monochrome camera is 

C = max/(X;U) = max[Ff(X) - H{X\Y)\ = maxh/(X) = Hc{X), (29) 

p{Xi) p{xi) p{xi) 

where I{X]Y) is the mutual information, Hc{X) = max if (X) the maximum information 

p{^i) 

entropy of source X for conventional imaging instrument, which is the Shannon Limit of 
the imaging system. According to the principle of maximum entropy^, the information 
content of an image is maximized when p{xi) is Gaussian distribution with average power 
constraints, which doesnt contain any useful information. Therefore, the entropy of the 
image with structured information H[X) has 

H{X) < Hc{X) = C. (30) 
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FIG. 7. The reconstructed spectral images of Mario &: Luigi with 30% sampling rate of 3D date- 
cube, showing all the channels from 520 ~ 580nm and 620 ~ 680nm. 


Eq. fl30|) shows that the image information acquisition efficiency of such conventional point- 
to-point imaging mode is lower than the Shannon Limit determined by Information Theory 
in optical imaging instruments. The channel capacity of an imaging system based on Infor¬ 
mation Theory for conventional optical imaging instruments is^^- 


C = Ndof log2(l + m), 


(31) 


where m is signal to noise ratio (SNR), Ndof is degrees of freedom and has 


Ndof = N, • N. ■ N. ■ N 


f • 1\ g • 1\ Q ' 


(32) 


where Nt,Ns, Nc and N^ are respectively time, spatial, color and polarization degrees of 
freedom. Spatial degrees of freedom Ng has^ 


,,, 0.61A,,, n.eiA, 1 

Ns = SW = {N^—){Ny—)- 


-fraclX/{2ay) = 1.22‘^N^Ny, 


(33) 


CXx Oiy A/ {2q.x) 

where S is the image area, W is the space bandwidth, ay and Nx, Ny are respectively 
the image-space aperture angle and the resolved pixel number in the image-space of coordi¬ 
nate X and y. The color degrees of freedom N^ depend on the number of spectral channels, 
while polarization degrees of freedom N^ is determined by the independent polarization 
state. According to Eqs. fl28|)fl^ and fl33ll . the channel capacity of the conventional camera 
in our experiment (where Nt^ = l,Nc^ = = l,Nx^ x Ny^ = 140 x 140, mi = 255) 

is Cl = A’dofi log 2 (l + mi) 2.33 x 10^, and the corresponding transmitted informa¬ 
tion of Fig,6(b) is Hi{X) ^ 1.07 x 10^ < Gi. In order to transmit the 520 


580nm 


wavelength ranges data, the required channel capacity of the conventional camera (where 
Nt^ = 1 ,Nc2 = 7,N^^ = 1,Nx2 X Ny^ = 140 x 140, m 2 = 255) is G ~ 1.63 x 10®, while 
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the required channel capacity in CISC spectral camera with 30% sampling rate in our ex¬ 
periment is C 3 ~ 4.90 X 10®. Ca < C 2 shows that CISC spectral camera has the higher 
information acquisition efficiency in a single exposure compared to the conventional camera. 
With the development of optical imaging technology, many new imaging technologies (such 
as CT imaged) are not based on the point-to-point imaging mode. However, because the 
correlation between pixels of image data doesnt be applied in the imaging reconstruction al¬ 
gorithm, the information acquisition efficiency of those new coding imaging technology also 
cant approaching the Shannon Limit determined by Information Theory for conventional 
optical imaging instruments. However, CISC imaging solution applies a spatial random 
phase modulation to satisfy the restricted isometry property (RIP)^ required by applying 
CS that makes the improvement of information acquisition efficiency of the imaging sys¬ 
tem possible. Comparing with CS imaging technology (such as Single-Pixel Imaging via 
Compressive Sampling^, coded aperture snapshot spectral imagers^), which forces on the 
compressive sampling of electric signal after photoelectric conversion to improve the chan¬ 
nel capacity utilization efficiency of the electric signal, CISC imaging solution improves the 
optical channel capacity utilization efficiency and achieves the compressive sampling of the 
image data during the imaging acquisition process, which opens the way of approaching the 
Shannon Limit determined by Information Theory in optical imaging instruments. As a new 
optical imaging technology, CISC spectral camera provides a unique solution for the spectral 
imaging of dynamic processes. This CISC imaging solution may also be expanded to other 
multi-dimensional information (such as polarization information) acquisition^, ultra-fast 
measurement^, and super-resolution imaging^'^'^. 


REFERENCES 

^E. Jacobs, Y. Fisher, and R. Boss, “Image compression: A study of the iterated transform 
method,” Signal Processing 29 , 251-263 (1992). 

^C. SHANNON, “A mathematical theory of communication,” The Bell System Technical 
Journal 27 , 397-423,623-656 (1948). 

®T. M. Cover and J. A. Thomas, Elements of information theory (John Wiley & Sons, 
2012) pp. 657-687. 

^P. Elias, “Optics and communication theory,” JOSA 43, 229-232 (1953). 


13 



®G. Francia, “Resolving power and information,” JOSA 45, 497-499 (1955). 

®W. Tan and Z. Wang, “Series of optical information theory,” Changchun Institute of 
Optics, Fine Mechanics and Physics, Chinese Academy of Sciences 1 (1964). 

^C. T. Di Francia, “Degrees of freedom of an image,” JOSA 59, 799-804 (1969). 

®W. Tan, “Optical information theory-retrospect and prospect,” Optics and Precision En¬ 
gineering 3, 17-22 (1982). 

®F. O. Huck, C. L. Tales, R. Alter-Gartenberg, S. K. Park, and Z.-u. Rahman, “Information- 
theoretic assessment of sampled imaging systems,” Optical Engineering 38, 742-762 
(1999). 

A. Strange, A. Duggins, W. Penny, R. J. Dolan, and K. J. Friston, “Information theory, 
novelty and hippocampal responses: unpredicted or unpredictable?” Neural Networks 18, 
225-230 (2005). 

I. Kolobov, Quantum imaging, Vol. 12 (Springer Science & Business Media, 2007) pp. 
79-110. 

H. Shapiro and R. W. Boyd, “The physics of ghost imaging,” Quantum Information 
Processing 11, 949-993 (2012). 

^^R. C. Gonzalez, R. E. Woods, and S. L. Eddins, Digital image processing using MATLAB 
(Pearson Education India, 2004) pp. 282-315. 

^^D. L. Donoho, “Compressed sensing,” Information Theory, IEEE Transactions on 52, 
1289-1306 (2006). 

^®E. J. Candes, J. Romberg, and T. Tao, “Robust uncertainty principles: Exact signal 
reconstruction from highly incomplete frequency information,” Information Theory, IEEE 
Transactions on 52, 489-509 (2006). 

C. Eldar and G. Kutyniok, Compressed sensing: theory and applications (Cambridge 
University Press, 2012). 

Wu, X. Shen, H. Yu, Z. Chen, Z. Liu, S. Tan, and S. Han, “Snapshot compressive 
imaging by phase modulation,” Acta Optical Sinica , 113-120 (2014). 

^®W. Gong and S. Han, “Super-resolution far-held ghost imaging via compressive sampling,” 
arXiv preprint arXiv:0911.4750 (2009). 

^®W. Gong and S. Han, “Experimental investigation of the quality of lensless super-resolution 
ghost imaging via sparsity constraints,” Physics Letters A 376, 1519-1522 (2012). 


14 



Wang, S. Han, and M. I. Kolobov, “Qnantnm limits of snper-resolntion of optical 
sparse objects via sparsity constraint,” Optics express 20, 23235-23252 (2012). 

Gong and S. Han, “High-resolntion far-field ghost imaging via sparsity constraint,” 
Scientific reports 5 (2015). 

Snn, M. P. Edgar, R. Bowman, L. E. Vittert, S. Welsh, A. Bowman, and M. Padgett, 
“3d compntational imaging with single-pixel detectors,” Science 340, 844-847 (2013). 

Zhao, W. Gong, M. Ghen, E. Li, H. Wang, W. Xn, and S. Han, “Ghost imaging lidar 
via sparsity constraints,” Applied Physics Letters 101, 141123 (2012). 

^^W. Gong, G. Zhao, J. Jiao, E. Li, M. Ghen, H. Wang, W. Xn, and S. Han, “Three- 
dimensional ghost imaging ladar,” arXiv preprint arXiv:1301.5767 (2013). 

^®W. Gong and S. Han, “Gorrelated imaging in scattering media,” Optics letters 36, 394-396 

( 2011 ). 

Bina, D. Magatti, M. Molteni, A. Gatti, L. Lngiato, and F. Ferri, “Backscattering 
differential ghost imaging in tnrbid media,” Physical review letters 110, 083901 (2013). 

^^O. S. Magaha-Loaiza, G. A. Howland, M. Malik, J. G. Howell, and R. W. Boyd, “Gom- 
pressive object tracking nsing entangled photons,” Applied Physics Letters 102, 231104 
(2013). 

^®W. Ghen and X. Ghen, “Object anthentication in compntational ghost imaging with the 
realizations less than 5% of nyquist limit,” Optics letters 38, 546-548 (2013). 

^®X. Xn, E. Li, H. Yu, W. Gong, and S. Han, “Morphology separation in ghost imaging via 
sparsity constraint,” Optics express 22, 14375-14381 (2014). 

Gheng and S. Han, “Incoherent coincidence imaging and its applicability in x-ray diffrac¬ 
tion,” Physical review letters 92, 093903 (2004). 

^^M. Zhang, Q. Wei, X. Shen, Y. Liu, H. Liu, J. Gheng, and S. Han, “Lensless fourier- 
transform ghost imaging with classical incoherent light,” Physical Review A 75, 021803 
(2007). 

^^H. Wang and S. Han, “Goherent ghost imaging based on sparsity constraint without phase- 
sensitive detection,” EPL (Europhysics Letters) 98, 24003 (2012). 

Zhang, Y.-H. Zhai, L.-A. Wu, and X.-H. Ghen, “Gorrelated two-photon imaging with 
true thermal light,” Optics letters 30, 2354-2356 (2005). 

^"^M. D’Angelo and Y. Shih, “Quantum imaging,” Laser Physics Letters 2, 567-596 (2005). 


15 



35X.-F. Liu, X.-H. Chen, X.-R. Yao, W.-K. Yu, G.-J. Zhai, and L.-A. Wu, “Lensless ghost 
imaging with sunlight,” Optics letters 39, 2314-2317 (2014). 

Giglio, M. Carpineti, and A. Vailati, “Space intensity correlations in the near held 
of the scattered light: a direct measurement of the density correlation function g (r),” 
Physical review letters 85, 1416 (2000). 

Cerbino, L. Peverini, M. Potenza, A. Robert, P. Bosecke, and M. Giglio, “X-ray¬ 
scattering information obtained from near-held speckle,” Nature Physics 4, 238-243 (2008). 

W. Goodman, Introduction to Fourier optics (Roberts and Company Publishers, 2005) 
pp. 154-160. 

^®C.-F. Cheng, D.-P. Qi, D.-L. Liu, and S.-Y. Teng, “The computational simulations of the 
gaussian correlation random surface and its light-scattering speckle held and the analysis 
of the intensity probability density,” Acta Physica Sinica 48, 1635-1643 (1999). 

^°A. Gatti, E. Brambilla, M. Bache, and L. A. Lugiato, “Ghost imaging with thermal 
light: comparing entanglement and classicalcorrelation,” Physical review letters 93, 093602 
(2004). 

W. G. Ghan, M. N. O’Sullivan, and R. W. Boyd, “High-order thermal ghost imaging,” 
Optics letters 34, 3343-3345 (2009). 

Gong and S. Han, “A method to improve the visibility of ghost images obtained by 
thermal light,” Physics Letters A 374, 1005-1008 (2010). 

^^F. Ferri, D. Magatti, L. Lugiato, and A. Gatti, “Diherential ghost imaging,” Physical 
review letters 104, 253603 (2010). 

Oymak, A. Jalali, M. Fazel, Y. G. Eldar, and B. Hassibi, “Simultaneously structured 
models with application to sparse and low-rank matrices,” arXiv preprint arXiv:1212.3753 
( 2012 ). 

Golbabaee and P. Vandergheynst, “Gompressed sensing of simultaneous low-rank and 
joint-sparse matrices,” arXiv preprint arXiv:1211.5058 (2012). 

Golbabaee and P. Vandergheynst, “Joint trace/tv norm minimization: A new efficient 
approach for spectral compressive imaging,” in Image Processing (ICIP), 2012 19th IEEE 
International Conference on (IEEE, 2012) pp. 933-936. 

^^H. Zhang, W. He, L. Zhang, H. Shen, and Q. Yuan, “Hyperspectral image restoration 
using low-rank matrix recovery,” Geoscience and Remote Sensing, IEEE Transactions on 
52, 4729-4743 (2014). 


16 



Hsieh, “Computed tomography: principles, design, artifacts, and recent advances,” 
(SPIE Bellingham, WA, 2009). 

F. Duarte, M. A. Davenport, D. Takhar, J. N. Laska, T. Sun, K. E. Kelly, R. G. 
Baraniuk, et ai, “Single-pixel imaging via compressive sampling,” IEEE Signal Processing 
Magazine 25, 83 (2008). 

^°D. Kittle, K. Choi, A. Wagadarikar, and D. J. Brady, “Multiframe image estimation for 
coded aperture snapshot spectral imagers,” Applied Optics 49, 6824-6833 (2010). 

^^S. P. Morgan and I. Stockford, “Surface-reflection elimination in polarization imaging of 
superhcial tissue,” Optics Letters 28, 114-116 (2003). 

Gao, J. Liang, C. Li, and L. V. Wang, “Single-shot compressed ultrafast photography 
at one hundred billion frames per second,” Nature 516, 74-77 (2014). 

^^D. L. Donoho, “Superresolution via sparsity constraints,” SIAM Journal on Mathematical 
Analysis 23, 1309-1331 (1992). 

^^E. J. Candes and C. Fernandez-Granda, “Towards a mathematical theory of super- 
resolntion,” Gommunications on Pure and Applied Mathematics 67, 906-956 (2014). 


17 



